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Abstract 

Sufficient conditions for synchronization of coupled Lienard-type oscil- 
lators are investigated via averaging technique. Coupling considered here 
is pairwise, unidirectional, and described by a nonlinear function (whose 
graph resides in the first and third quadrants) of some projection of the 
relative distance (between the states of the pair being coupled) vector. 
Under the assumption that the interconnection topology defines a con- 
nected graph, it is shown that the solutions of oscillators can be made 
converge arbitrarily close to each other, while let initially be arbitrarily 
far apart, provided that the frequency of oscillations is large enough and 
the initial phases of oscillators all lie in an open semicircle. It is also 
shown that (almost) synchronized oscillations always take place at some 
fixed magnitude independent of the initial conditions. Similar results are 
generated for nonlinearly-coupled harmonic oscillators. 

1 Introduction 

Synchronization in coupled dynamical systems has been a common ground of 
investigation for researchers from different disciplines. Most of the woA in the 
area studies the case where the coupling between individual systems is linear; 
see, for instance, [HJ [TH [T3J [501 13 HH]- Nonlinear coupling is also of interest 
since certain phenomena cannot be properly modelled by linear coupling. A 
particular system exemplifying nonlinear coupling that attracted much attention 
is Kuramoto model and its like [51 |T3] . Among more general results allowing 
nonlinear coupling are [21 I17j where passivity theory is employed to obtain 
sufficient conditions for synchronization under certain symmetry or balancedness 
assumptions on the coupling graph. 

In this paper we study the synchronization behavior of an array of nonlinear 
planar oscillators. We let the individual oscillators share identical dynamics and 
the coupling between them be nonlinear. We consider Lienard-type oscillators 
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that have been much studied due to their close relation to real-life systems and 
applications [TJ |6l [16] . A particular example is van der Pol oscillator [71 [11] . 

What we investigate here is the relation between frequency of oscillations and 
synchronization of the oscillators forming the array. The array is formed such 
that some of the oscillators are coupled to some others via a nonlinear function. 
Coupling considered is of partial-state nature. That is, if an oscillator affects 
the dynamics of another, the associated coupling term is not a function of the 
state vector of the oscillator that is affecting, but only of some projection of that 
vector. We make no symmetry nor balancedness assumption on the coupling 
graph. 

Our finding in the paper is roughly that (almost) synchronization occurs 
among the oscillators (at some magnitude independent of the initial conditions, 
the coupling, and the frequency of oscillations uj) if the following conditions 
hold: (a) the frequency of oscillations is high, (b) there is at least one oscil- 
lator that directly or indirectly affects all others, and (c) the initial phases of 
oscillators, when each is represented by a point on the unit circle, all lie in an 
open semicircle. More formally, what we show is that if the coupling graph is 
connected and the initial phases of oscillators lie in an open semicircle, then 
the solutions of oscillators can be made converge arbitrarily close to each other, 
while initially being arbitrarily far from one another, by choosing large enough 
oj. Incidentally, as sort of a byproduct of our analysis for nonlinear oscillators, 
we also generate a similar result for an array of nonlinearly-coupled harmonic 
oscillators. We show that harmonic oscillators (almost) synchronize provided 
that the coupling graph is connected and the frequency of oscillations is high. 
Different from nonlinear oscillators, the initial phases of harmonic oscillators 
do not have any effect on synchronization, at least when uj is sufficiently high. 
Also, again unlike nonlinear oscillators, synchronized oscillations can take place 
at any magnitude, depending on the initial conditions. 

In establishing our results we use tools from averaging theory |3] . Our three- 
step approach is as follows. We first apply a time- varying change of coordinates 
to the array, which keeps the relative distances between the states of oscillators 
intact (Section [1]). This change of coordinates renders the array periodically 
time-varying. Then we take the time-average of this new array and show, for 
the average array, that the oscillators synchronize if their initial phases lie in 
an open semicircle and the coupling graph is connected (Section [5]) . Finally we 
work out what that result means for the original array. Namely, we show that 
the higher the frequency of oscillations the more the original array behaves like 
its average (Section [B]). 

2 Preliminaries 

Let K>o denote the set of nonnegative real numbers. Let | • | denote Euclidean 
norm. A function a : M>o — >■ K>o is said to belong to class-/C {a £ /C) if it is 
continuous, zero at zero, and strictly increasing. It is said to belong to class-/Ctx> 
if it is also unbounded. Given a closed set S C M" and a point x G K", |a;|5 
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denotes the (Euclidean) distance from x to S. Convex hull of S is denoted by 
coiS. Number of elements of a finite set X is denoted by #1. 

A [directed) graph is a pair (N, £) where A/" is a nonempty finite set (of nodes) 
and £^ is a finite collection of ordered pairs (edges) (n^, rij) with n.;, iij e TV. 
A directed path from ni to n; is a sequence of nodes (ni, . . . , ni) such that 
{ni, rii+i) is an edge for iG{l, — 1}. A graph is connected if it has a node 

to which there exists a directed path from every other node0 

A set of functions {7^^ : M — ?► R}, where i, j — 1, . . . , m with i ^ j, describes 
(is) an interconnection if the following hold for all i, j and all s G R: 

(Gl) 7,y (0) = ands7,y (s) >0. 

(G2) Either %j{s) = or there exists a ^ K, such that |7ij(s)| > a(|s|)- 

To mean 7y (s) = we write 7^ = 0. Otherwise we write 7^^ ^ 0. The graph 
of interconnection {7,., } is pair (TV, £), where TV = {ni, . . . , n™} and £ is such 
that [ui, Uj) € £ iff 7^ ^ 0. An interconnection is said to be connected when 
its graph is connected. To give an example, consider some set of functions 
G := {7ii : h j = 1, ■ • • , 4}. Let 713, 723, 724, 732 be as in Fig. 1 and the 
remaining functions be zero. Note that each 7^^ satisfy conditions (Gl) and 
(G2). Therefore set Q describes an interconnection. To determine whether Q is 
connected or not we examine its graph, see Fig. 2. Since there exists a path to 
node 714 from every other node, we deduce that the graph (hence interconnection 
Q) is connected. 

Consider two half lines on (originating from the origin) such that the 
(smaller) angle between them is strictly less than tt. Then their convex hull C is 
called a cone, see Fig. 3. Angle of cone C, denoted by ZC, is the angle between 
the two half lines. Note that ZC 6 [0, tt) and a half line is a cone with zero 
angle. Given r > 0, a set Wr C is called an r-wedge if it can be written as 
Wr — co(C n {a; G R^ : ri < |a;| < r2}) for some cone C and < ri < r < r2- 
See Fig. 4. 

3 Problem statement 

A general model for a planar oscillator is given by Lienard's equation 

q + f{q)q + g{q)=0. (1) 

Sufficient conditions have been established on functions / and g in order for 
system ([T]) to admit a unique, stable limit cycle encircling the origin of the 
phase plane; see, for instance, [T] Thm. 1]. In this paper we adopt those gen- 
eral conditions on function /, which are stated later in the section. Regarding 
function g, we assume linearity. Namely, we study the case g{q) — uj'^q, where 
uj is some constant. Although this assumption is restrictive compared to what 
is assumed in Thm. 1], it nevertheless still allows us cover some physically 

^This is another way of saying that the graph contains a directed spanning tree. 
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7l3(s) 723(s) 




Figure 1: Some examples of coupling functions. 



important cases. Most famous example that falls into the class of systems being 
studied here would be van der Pol oscillator [8j. In this paper we search for 
sufficient conditions that yield synchronization of a number of coupled Lienard 
oscillators. Below we give the precise description of the problem. 
Consider the following array of coupled Lienard oscillators 



f{qi)Pt 



(2a) 
(2b) 



where w > is the frequency of oscillations and {7^^} is a connected inter- 
connection. Let G denote the state of «th oscillator, i.e., £,i = [qi Pi]^ ■ 
Sometimes we choose to handle this array ([2|) of to planar oscillators as a single 



system in . 



If we let e . . . C 



TiT 



and 



ujpi 



-^qi - fiqi)pi + E7ii(pj -Pi) 



LOp„ 



-Uqm - fiqm)Pm + J2 [Pj - Pm) 
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Figure 2: Graph of interconnection Q. 




Figure 3: A cone C. 



then array ^ defines the below system 

i = f-iC, u;) . (3) 

We assume throughout the paper that jij and / : R — !> R are locally Lipschitz. 
Letting F{s) J^^ f{(T)da, we further assume the following. 

(LI) / is an even function. 

(L2) F{so) = for some so > 0; F is negative on (0, sq); P is positive, nonde- 
creasing, and unbounded on (sq, oo). 

In this paper the question we ask ourselves is the following. 

What is the effect of lo on the synchronization behavior of array ^ ? 

Our approach to the problem is as follows. We first apply a norm-preserving, 
time-varying change of coordinates to array ([2]) and obtain a new array whose 
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Figure 4: An r- wedge W^. 



righthand side is periodic in time. Since this coordinate change does not aher 
the relative distances between the states of oscillators, it is safe to look at the 
new array to understand the synchronization behavior of the original one. Hence 
we focus on the new array. Although the righthand side of the new array does 
not look any pleasanter than the original one, it nevertheless has an average 
since it is periodic in time. Once this average is computed we realize that it is 
very simple to understand from synchronization point of view. From then on we 
start going backwards. Using averaging theory, we establish what our findings 
on the average array imply, first for the time-varying array and eventually for 
the original array. 



4 Change of coordinates 



We define 8(10) e K^^^ and H,V eR^ 



as 



io 
uj 



:= [1 0] , V ~ [0 1] . 



Then we rewrite ([2]) as 



(4) 
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We can by write 



Let x,{t) := e-'S'(")*f,(i) and x := [xf 



= —f {[cos Lut sin ujt\xi) 



- sin cji 
coscji 



E 



— sin ojt 
cos wt 



— sincji cosujt\xi 
7ij([— sincji cosa;i](a;j — x^)) 



(5) 



Let 



i'^'^ix, UJt) := 



el'^ix, cut) 



Then system in M^™ corresponding to array (O reads 



(6) 



Remark 1 Since the change of coordinates is realized via rotation matrix e^"^^")*, 
the relative distances are preserved, that is \xi{t) — Xj(t)\ = \S,i{t) —^j(t)\ for all 
t and all i, j . This means from synchronization point of view that the behavior 
of array ([2]) will be inherited by array ([5|) . 

Exact analysis of ([5]) seems far from yielding. Therefore we attempt to under- 
stand this system via its approximation. 

5 Average array 



tions / : — ^ and jij : — > are given by 



Observe that the righthand side of ([5]) is periodic in time. Time-average func- 

sin</7 cos ip\vdip 



2-K 



fiv) ■= ^ f {[cos if smip]v) 



- simp 
cos if 



and 





— simp 


Jo 


COS(^ 



jij{[—simp cosip]v)dip. 



Then the average array dynamics read 



Vi = -fiVi) + E '^v^^J ~ • 



(7) 



(8) 
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Let 77 := [yyf . . . 77 Jj^ and 

-fim) + J2lij{V3 -Vi) 

-fiVrn) + E ImjiVj ^ Vni) _ 

Then system in M^™ corresponding to array ([8]) reads 

,) = F«(,y). (9) 
Remark 2 A^oie i/iaf instead of array ^ if we start with the below array 

where C G M^^^ satisfies CC^ = 1, we still reach the same average array ([5]). 
Therefore for the analysis to follow the coupling need not be through velocities 
as is the case in ([2]). Any projection of the state CS,i is as good as any other for 
coupling. For instance, the results in this paper hold true for the below array 

(ji = ojp, + ^'yij{qj ~ Qi) 

Pi = -cj^fj - f{qi)Pi ■ 

Theory of perturbations |3l Ch. 4 § 17] tells us that, starting from close initial 
conditions, the solution of a system with a periodic righthand side and the solu- 
tion of the time-average approximate system stay close for a long time provided 
that the period is small enough. Therefore ^ should tell us a great deal about 
the behavior of ([5|) when cj 3> 1. Understanding ([8|) requires understanding 
average functions / and jij. Let us begin with the former. 

V 

Lemma 1 We have f{y) — r where k : M>o — >■ M is 



<^)--=l [ fi^)\l^-^da. (10) 
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Proof. Given v G M^, let r = \v\ and 9 G [0, 27r) be such that r [cos 61 sm^]-^ = v. 
Then, by using standard trigonometric identities, 



1 

2^ 



/(r(cos(/3cos0 + sin sin 6*)) 



X r (— sin ip cos 6* + cos cp sin ( 



- sm If 
cos 



27r 

= ^TT I f{r{cos{Lp~e))sm{(p-0) 

2,71 



r 
r 



f{r cos if) sin 



- sm{(p + 6) 
cos{ip + 6) 



- sm ip 

cos if 

dip 



dip 



2tt 



— I f{r COS p)amp 



r 

2^ 



/(r cos p) BIT? pdp 



sm p cos u — cos (/3 sm ( 
cos p cos 6 — sin (/3 sin ( 

cos 6 



dp 



sa\9 



^ / f {r comp) simp cos pdp 



smt/ 
cos 9 



where the second term is zero since f{r cosp) sin 293 is an odd function and we 
can write 



/(r cos p) sin ip cos ipdp 



f{r cos ip) sin 2pdip . 



Therefore 



27r 



/(w) = ( ^ / f{rcosp)s\TL ipdp 



cos 9 
s\\\9 



/ /(rcos(p)rsin2v"^(^ ) . (11) 
27r7_^ / \v\ 



Since f (r cos p)r sva^ ip and /(r sin 93)r cos^ p are even functions we can write 



1 

2^ 



/(r cos iy9)r sin^ pdp 



f{r cos (p) r sin pdp 

Tr/2 



7r/2 



/(r sin (/3)r COS (/3d(/5 



2 /-"/^ 



/(r sin iy9)r cos^ pdp . 



Then by change of variables a := rsinp we obtain 



1 

2^ 



f{r cos p)r sin^ 1^9 dp 



f{a)^l-—da. 



(12) 



Combining dTTI) and gives the result. 
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Claim 1 Defined in (|10p . function k satisfies (L2). 

Proof. Recall that F{s) — f satisfies (L2). Therefore there exists Sg > 
such that F{so) = 0, F{s) < for s G (0, sq), F is positive, nondecreasing 
on (sq, oo), and F{s) -> oo as s — oo. Given s > 0, let Sg ■ [0, s] [0, 1] 
be 6s{(y) :— \/l — (cr/s)^. Then we can write k{s) — ^ f{a)Ss{(7)da. We 
observe the following properties. 

(Dl) Given s > 0, function 6s{-) is strictly decreasing. 

(D2) Given s > t > 0, map a M> Ss{cr) — St{(j) is strictly increasing on [0, t]. 
(D3) For each t > there exists s > t such that 5s{(j) > i for u G (0, t). 
It follows from integration by parts that 

/ fia)S,{a)da <0 Vs G (0, so) • (13) 
Jo 

Let us convince ourselves that ([T3| is true. Let Sg{a) :— dSs{<j)/da. Given 
s G (0, So) can write 



f{a)S,{a)da = F{a)6,{a) ~ / F{a)5',{a)da 

Jo 

F{a)Ua)da 



< 

since for a G (0, s) we have (5;(cr) < by (Dl) and F{a) < by (L2). 
Next we show 

/•So psa 

/ f{a)6s{a)da> / f{a)St{a)da for s > t > Sq . (14) 

Given s > t > sq, let A((t) := (5s((t) — ^t(c) and A'((t) := dA{a)/d(T. We can 
write 



/(a)^.(a)d(7 = / f{a)A{a)da+ / /(a)5t(a)da 

Jo JQ 

= Fia)A{a)'° - f " Fia)A'{a)da+ I f{a)5t{<j)da 



F{a)A'{a)da+ [ ° f{a)St{<j)da 



> / fia)5t{a)da 
Jo 

since for a G (0, so) we have A'(ct) > by (D2) and F{a) < by (L2). 
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Since F is positive, nondecreasing on (sq, oo), we have /(z) > for z > sq. 
Therefore (D2) implies 



f{a)6s{a)da > / f{a)6t{(y)da for s > t > sq . 



(15) 

Combining (1141) and ([15]) we obtain 

f{(T)Ss{<j)da > / f{a)5t{(T)da for s > i > sq . (16) 
Jo 

Moreover, since -F(s) — oo as s — > oo, (D3) readily yields 

lim / f {a)5 s{(j)da = oo. (17) 
Combining dTS]), dH]), and ([171) gives the result. ■ 
Lemma 2 We /lawe ^ij{v) = PijUvDy-^ where pij : R>o — > K>o *s 

1 f^"" 

P^ji^) TT Jij{ssiinp)siinpdip . (18) 

Proof. Given w e M^, let r = |w| and 9 € [0, 27r) be such that r[— cos 61 sin0]^ = 
i>. Then, by using standard trigonometric identities, 



2tt 







1 

1 

2^ 
1 

2^ Jo 
-f 



- sm(p 
cos (p 

- sa\ip 
cos ip 



— sin((p — 6*) 
cos((p — 6*) 



7ij (r(sin ip cos + cos p sin 9))dip 
7ij(rsin((y9 + 6'))d(/5 
7i-,(r sin93)d(^ 



sm iy9 cos U + cos sm ( 
cos if cos + sin tp sin ( 



7ij (r sin 93) sin pdp 



2-iT 



jij{r sin (^)d(p 



cos 6* 
sin6' 



~^ \2tt J 7ii(''sin(/j)cos(y9(i<^ 
We focus on the second term in ([T9| . Observe that 



smt/ 
cos 



(19) 



7y (r sin 1^9) cos pdp 



-7r/2 



7»j 



(^r sin + I) ) cos ((^ + |) dp 
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since the integrand is an odd function on the interval of integration. Likewise, 



27r 



7ij (r sin 93) cos ipdip 



-it/2 



-Tr/2 



37r\ 

7ij ( r sm ( + — 1 ) cos ( (/3 + 



Therefore 

/.27r 

7y (r sin (/s) cos 

7y- (r sin (/j) cos ipdip - 

I 

= 0. 

Combining (IT51) and ([^0]) we obtain 

^27r 



7ij (r sin (p) cos ipdp 



(20) 



7ij(i') = 



1 
2^ 



7ij (r sin tfi) sin 1^9^93 



— COSf 

sin( 



Hence the result. 



Claim 2 Defined in (|18p . Junction Pij{s) = iff ^ij — 0. Otherwise, there 
exists a G IC such that Pij{s) > a{s) for all s. 

Proof. That pij{s) = if jij = is evident from the definition. Suppose 
7ij ^ 0. Then by (G2) there exists ai G /C such that |7ij(s)| > ai(|s|). Then 
we can write 



Py(s) 



1 /-^^ 

— / '-f.ij{ssmp)smpdp 
Jo 



> 



1 

2^ 
1 

2^ 



|7ij (s sin (^) 1 1 sin p\dp 



ai{s\ siniy9|)| sinp\dp =: 02(3) . 
Note that a2 is a class-/C function. Hence the result. 



In the remainder of the section we generate two results on average array (|8|. 
In the first of those results (Theorem[T|) we establish that each of the oscillatort|j 
in ([8]) eventually oscillates with some magnitude no greater than some constant 
p which only depends on function / and not on initial conditions. In the sec- 
ond result (Theorem [5]) we assert that if the initial conditions are right then 



■^We admit that the word oscillator may have been an unfortunate choice here since there 
is no oscillation (in the standard meaning of the word) taking place in the average array. 
Averaging gets rid of oscillations. 
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the oscillators eventually synchronize both in phase and magnitude, where the 
magnitude equals p. Initial conditions' being right roughly corresponds to the 
following condition. If we depict each oscillator's initial phase with a point on 
the unit circle then those points should all lie in an oper|f| semicircle. Constant 
p we mentioned above is indeed defined as the unique positive number satisfying 

K{p) = . 

Existence and uniqueness of p is guaranteed by Claim [TJ Note then that the 
first result is about the asymptotic behavior of the solutions of system ^ with 
respect to the following set 

f^:={r;eR2™ : |r;,| <p}. 

Likewise, the second result has to do with the asymptotic behavior of the solu- 
tions of system ^ with respect to 

7^ {r; e ]R2m . ^ \t],\=p for all i, 3) . 

Theorem 1 Consider system Set B is globally asymptotically stable. 

Proof. We will establish the result by constructing a Lyapunov function. 
Claim [1] implies that there exists ai & K, such that ai{s ~ p) < k{s) for s > p. 
Then we can find a2 £ /C such that 

a2(|?7|B) < ai (^max{0, max{|?7i| - p}}j max|?7i| . 
Let our candidate Lyapunov function V : M^'" — > R>o be 

V{r]) i max{0, max{|?7,p - p^}} . 
Then there exist as, 0:4 G /Coo such that 

a3{\v\B)<V{fi)<a4{\v\B)- (21) 

Now, given some 77 with |?7|e > 0, let (nonempty) set of indices I be such that 
^(hiP ~ P^) — ^iv) iff * £ ^- Observe that for i 6 I, point r]i lies on the 
boundary of the smallest disk (centered at the origin) that contains all points 
r]j . Therefore i € X implies rjf {rjj — rji) < for all j . Then by Lemma [T] and 

•^That is, cndpoints are not included. 
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Lemma [5] we can write (almost everywhere) 



< 

< 
< 

Result follows from ([HI) and ^^j. ■ 

Remark 3 Theorem{J\ does not require interconnection {"fij} to be connected. 
This is clear with (or even without) the proof. 

Theorem 2 Consider system Q . Set TZ is locally asymptotically stable. In 
particular, co{77i(0), ?7m(0)} n {0} = implies r/{t) -^77* as t (x for 
some r/* G TZ. 

To prove the theorem we use an invariance principle (similar to that of 
LaSalle's) for which we need to tailor an invariant set for our system. To this 
end we introduce some notation. Given 77 e R^™ with co{?7i, . . . , r/m} n {0} = 
we let C{ri) C and Wp(r/) C respectively be the smallest cone and smallest 
p-wedge containing set {771, . . . , ri„i}. (Recall that p > satisfies k{p) — 0.) 
Fig. 5 depicts Wp{ri). Note that 

Wpiv) C Civ) . (23) 

For C e R^™, when we write C e ^^(77)" we mean that Q G Wp (77) for all 
i = 1, . . . , m. Note that 

C e >Vp(r7)™ =^ >Vp(C) C >Vp(77) . (24) 

Proof of Theorem [2l We first show local stability of TZ. Let 77 G K^™ be such 
that co{77i, . . . , 77,,,, } n {0} = 0. Consider the first term of the righthand side 
in (|5]). Lemma [1] and Claim [T] tell us that vector ~f{rji) (tail placed at point 
77i) points toward the origin if [77^1 > p and away from the origin if \rii\ < p. 
Therefore, for 77^ on the boundary of Wp(7;), vector —f{r]i) never points outside 
Wp (77). Now consider the sum term in (|8]). Lemma [2] and Claim [2] tell us that 
each summand — rji), if nonzero, is a vector pointing from rji to rjj. Hence 



max 77f I -/(77,) + ^7ii('7i " Vi) 



max -iiffivi) + V vfjijiVj ^ Vi) 
max I -K(|77.|)|7;..|+^ ^'f'^^~^'^ f (77,-770 



-K ^niax |77i I J max \rii \ 
-ai (^max{|77i| - r}^ max \r]t\ 

-M\v\b)- (22) 
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Figure 5: An example with 77 G R^°. Dots represent the positions of on plane. 
The shaded region is Wp (77). 



sum lij cannot be pointing outside convex hull 00(771, . . . , 77^} for 

77i on the boundary of the convex hull. Since 00(771, . . . , C Wp (77) we deduce 
therefore that the righthand side of dS]), that is, vector —f{rii)+J2j=/.i lij(.''lj~Vi)j 
never points outside Wp (77) for 77^ on the boundary of Wp(77). Hence compact 
set Wp(77)'" is forward invariant with respect to system ([9|). Observe also that 
for 77* G 7^ we have Wp(7y)'" {v*} as 77 — !> 7/*. Hence set TZ is locally stable. 
Next we show attractivity. 

Let G denote the graph of interconnection {7ij}. Note that G is connected 
by assumption. By Lemma [2] and Claim [2] we know that pij is continuous and 
zero at zero. Moreover, if there is no edge of G from node i to node j then 
pij{s) = 0. Otherwise there exists a € K. such that Pij{s) > a(s) for all s. 

Consider system dH). Let 00(771(0), . . . , 77^(0)} n {0} = 0. Due to ([24]) and 
forward invariance of Wp(77(t))'", for t2 > ti > we can write Wp{ri{t2)) C 
Wp(77(ti)). By (|23)) map t /LC{-q{t)) is hence nonincreasing. It is also 
bounded from below by definition. Therefore there exists 6 G [0, tt) such that 
limt_j>oo ^C{r]{t)) = 9. We now claim that 9 = 0. 

Suppose not. Then, by continuity, there must exist a solution to sys- 
tem ® with C(0) G Wp(7;(0))'" such that /C(C(t)) = 9 for all t > 0. Recall 
that Wp (C(0)) is forward invariant. By therefore 

C(C(i))=C(C(0)) Vt>0. (25) 

Let Hi and H2 be two half lines convex hull of which equals C(C(0)). Since 
61 > 0, we have Hi ^ H2. Then we observe that {Ci(0), . . . , Cm(0)} n 7^ 
for i = 1, 2. Let 77,1, . . . , rim denote the nodes of graph G. Since G is connected 
we can find a node ni to which there exists a path from every other node. 
Now, without loss of generality, assume that 0(0) ^ ^i- Let Ti{t) denote 
the set of indices at time < > such that Q{t) G Hi iff i G Xi{t). Clearly, 
I ^ 2^1(0) and 1 < #li(0) < 77i — 1. Hence connectedness of G implies that 
there exist iq G Ii(0) and jo G {!:•■•: "i} \ li(0) such that (71^,-,, 77,^^) is an 
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edge of G. Existence of edge {ni„, Uj^) implies that ^ 0. Now, observe 

that Cjo(O) - Cio{0) 7^ since Cio(O) G Tii and Cio(O) ^ "^i- % Lemma [3 
7ioio(Cjo(0) ~ Cio(O)) is a nonzero vector that is not parallel to Hi- By ([S|) we 
can write 

6o(o) = i^ + 7»ojo(Oo(o)-C.o(o)) 

where 

« = -/(C.o(0))+ E 7.o.-(0(0)-Go(0)). 

J^^o Jo 

By earlier arguments (see the first paragraph of the proof) we know that vector 
V (tail placed at point Ci(0)) cannot point outside C(^(0)). Therefore vector 
(■20(0) is nonzero and not parallel to Hi- This lets us be able to find > such 
that Cio(^i) ^ ^1 ^^'^ #2^1 (^i) < #2^1(0) — 1- We can continue the procedure 
and find a sequence of instants tfc > such that #Xi(tfe) < #li(0) — fc. Since 
#Ji(0) is finite, there should exist t* > such that #Ii(t*) = 0, that is, 
Ii(t*) = 0. This observation translates to "Hi n {Ci(i*), • • • , Cm(^*)} = 0, which 
yields C(C(i*)) 7^ C(C(0)) which contradicts with ([25]) . Hence our claim is valid 
and limt^oc ^C{rj{t)) = 0. 

That Z.C{-q{t)) means that C(77(i)) U for some half line U. We 
also know that ri{t) e >Vp(7?(0))" for alH > and, by Theorem [fl ry(t) B. 
Therefore solutions 'qi{t), for i = 1, . . . , m, converge to the following set 

5 := 7^ n Wp(?7(0)) n {u e : |w| < p} . 

Note that 5 is a line segment that satisfies S — {\u : 5 < \ < p} ioi some unit 
vector M G and 5 > Q. Solution 77(-) will converge to the largest invariant 
set in S"^ . (Note that 5™ itself is forward invariant.) Take any solution (,{■) to 
system © with C(0) G 5™. At any given time i > let j be such that point 
C,i{t) is farthest from the point pu. Then by Lemma [TJ Claim [1] Lemma [51 and 
Claim [5] we have 

/ if Ut)=PU 

'''^^> '{ Xu if Q{t)^pu 

for some A > —K{\C,i{t)\) > 0. Therefore max^ \Q{t) — pu\ Q as t ^ 00. We 
then deduce that {pit}™ is the largest invariant set in 5™. Hence the result. ■ 

Remark 4 In Theorem\^we require condition co{?7i(0), . . . , ?7m(0)} H {0} = 
for synchronization. If we express each rji in polar coordinates {ri, 6i), where 
fi > and 9i G [0, 27r), then the required condition is equivalent to that ri[Q) ^ 
and 9i{0) € {9* , 9* + tt) for all i and some 9* . In other words, if we denote 
each phase 9i{0) by a point on the unit circle, then the condition is equivalent 
to that they all lie in an open semicircle. Is this condition really essential for 
synchronization ? Theoretically speaking, yes. That is, one can easily find initial 
conditions that make an equilibrium but are not on the synchronization manifold. 
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However those equilibria may be practically disregardable if they are unstable, 
since slightest disturbance would rescue the system from being stuck at those 
points. Now we ask the second question. Does there exist a stable equilibrium 
outside synchronization manifold? This we find difficult to answer. We would 
nevertheless like to report that our attempts to hunt one via (limited) simulations 
on coupled van der Pol oscillators have failed. 

Let us recapitulate what has hitherto been done. We have begun by an array 
of coupled Lienard oscillators ^ — ^(^. a;) with frequency of oscillations uj. By 
looking at it from a rotating (at frequency a;) reference frame we have obtained 
new array x = £^''^{x, ujt) whose righthand side is periodic in time. Exploiting 
this periodicity, we have then computed average dynamics rj = i^^{ri). Finally 
we have shown for the average array that if oscillators' initial phases all reside in 
an open semicircle then they synchronize both in phase and magnitude, where 
the magnitude should equal p. 

In the next section, based on our findings on the synchronization behavior 
of average array, we show that oscillators of array x — £^'-^{x, tut) and therefore 
of array ^ — £{^, uj) get arbitrarily close to synchronization for lo large enough. 
We use tools from averaging theory to establish the results. 

6 Synchronization of Lienard oscillators 

The following definition is borrowed with slight modification from |18) . 

Definition 1 Consider system x ~ g{x, uj, t). Closed set S is said to be 
semiglobally practically asymptotically stable (with respect to u) if for each 
pair (A, 5) of positive numbers, there exists w* > such that for all lu >uj* the 
following hold. 

(a) For each r > 5 there exists £ > such that 

\xito)\s <e =^ \x{t)\s <r Vt > to • 

(b) For each e < A there exists r > such that 

\x{to)\s<£ =^ \xit)\s<r \ft>to. 

(c) For each r < A and e > S there exists T > such that 

\xito)\s<r =^ \xit)\s<e Vt>to+T. 

Result [m Thm. 2] tells that the origin of a time- varying system is semiglob- 
ally practically asymptotically stable provided that the origin of the average 
system is globally asymptotically stable. The analysis therein, with almost no 
extra effort, can be extended to cover the case where the attractor is not a 
singleton but only compact. Note that set B is compact. Theorem [T] therefore 
yields the following result. 
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Theorem 3 Consider system Set B is semiglobally practically asymptoti- 
cally stable. 

The following result is an adaptation of a general theorem on averaging \TE[ 
Thm. 2.8.1]. 

Lemma 3 Let map g : M" x M — > R" be locally Lipschitz and satisfy g{x, (f + 
2ir) — g(x, ip) for all x G M" and G R. Define average function g : M" — >■ K" 
as 

In 



1 r'' 

271- Jo 



Let x{-) and rj{-) denote, respectively, the solutions of systems x = g{x, cot) and 
fj = g{f]), where w > 0. Then for each compact set 2? C M" and pair of positive 
real numbers (T, e) there exists w* > such that if 

• id > LO* , 

• 77(0) — x{Q), and 

• x{t) e V for all t e [0, T] 

hold then \ri[t) — x{t)\ < e for all t e [0, T\. 

Lemma 13] lets us establish the following result. 

Theorem 4 Consider system (|6]). For each pair (V, S), where V C is a 
compact convex set that does not include the origin and S > 0, there exists 
UJ* > such that for all ui > u* the following hold. 

(a) There exists £ > such that 

k(0)|7j <e =^ \x{t)\n<5 > . 

(b) There exists r > such that 

x(0) e 2?" =^ \x{t)\n <r Vt > . 

(c) There exists T > such that 

a;(0)er'" => \x{t)\n<5 yt>T. 

Proof. Let us be given (V, 5). Let x{-) and ry(-) respectively denote the solu- 
tions of system (|6]) and system ([9]). We first work on part (a). Choose some 
61 > such that \x\tz < Si implies co{a;i, Xm} D {0} = 0. Then set 
82 ■— ^min{(5i, 6}. Choose e e (0, ^2) such that 

h(0)k<e ^ \v{t)\n<S2 Wt>0. 
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Such e exists since TZ is locally asymptotically stable for system ^ by Theo- 
rem[2l Moreover, |77(0)|tc < £ implies co{?7i(0), . . . , ''7m(0)}n{0} = 0. Therefore 
we can find Ti > such that 

\rimn<e =^ hWk < | Vt > Ti . 

Let ei := min{|, 52}- By Lemma [3] there exists wi > such that \x{t) — ri{t)\ < 
El for all t S [0, Ti] provided that uj > uji, \x{t)\Ti < 6 for all t G [0, Ti], and 
x{0) = 7/(0). Let now uj > coi and \x{0)\n < e. Set r]{0) = x{0). Then for all 
t & [0, Ti] we can write 

\xit)\n < \ii{t)\Tz + \x{t)-rjit)\ 

< S2 + ei 

< 5. 

Also 

\x{Ti)\n < \ii{T,)\n + \x{T,)~Tj{Ti)\ 
e 

< £• 

We can repeat this procedure for the following intervals [kTi, (fc + l)Ti], k = 
1, 2, . . . Therefore \x{0)\tz < e implies \x{t)\']z < 5 as long as w > oji. 

Part (b). By Theorem|3]we can find uj2 > and ri > such that uj > u!2 and 
x(0) G 2?™ imply < ri for all t > 0. Then we can find r > such that 

\x\b ^ "f"! implies \x\fi < r. Therefore uj > lu2 and a;(0) G 2?™ imply |x(t)|7j < r 
for all t > 0. 

Part (c). First recall, from the proof of Theorem HI that any p- wedge is 
forward invariant with respect to system ©. Choose a p- wedge Wp such that 
V C Wp. Note that Wp is compact, convex, and does not include the origin. 
By Theorem |4] there exists T > such that 

^(0)GW™ ^ h(r)k<£i 

where ei is as found in part (a). Then define S C M^™ as 

5 := {x + eiv:xe W™, |u| < 1} . 

By Lemma |3] there exists ui^ > such that \x{t) — ri{t)\ < ei for all t G [0, T] 
provided that w > wa, G S for all i G [0, T], and ^(O) = ?7(0). Let now 
oj > UJ3 and x(0) G P™. Set 7?(0) = a;(0). Then for all t G [0, T] we have 
k(^) ~ '^l*)! ^ £i- Therefore 

< \v{T)\n + \x{T)~f^{T)\ 

< 61+ ei 

< e. 
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From part (a) we know that \x{T)\ii < e implies \x{t)\-ji < S for all t >T, which 
was to be shown. 

We complete the proof by setting u* := max{a;i, W2, W3}. ■ 

We next present the below result for the array of coupled Lienard oscillators. 
This result directly follow from the observation mentioned in Remark [T] and the 
fact that sets B and TZ are invariant under rotations in M^. 

Theorem 5 Theorem and Theorem [7] hold true for system ^ . 

Theorem [3] roughly says that solutions Xi{-) of oscillators of array ([5|) can 
be made converge to an arbitrarily small neighborhood of the disk {u G : 

^ p} starting from arbitrarily large initial conditions by choosing oscillation 
frequency u arbitrarily large. Theorem |4] says that by choosing uj arbitrarily 
large, solutions Xi{-) can be made eventually become arbitrarily close to each 
other and to the ring {v € : \v\ = p}, starting from within an arbitrary 
convex compact set that does not contain the origin. Finally Theorem [5] says 
that these two results should hold also for array of Lienard oscillators ^ . 

Theorem [5] marks the end of our answer to the question that we asked for 
coupled Lienard oscillators in Section [3l It is hard not to realize that the basic 
idea forming the skeleton of this answer serves as a solution approach also for 
the problem of understanding the synchronization behavior of coupled harmonic 
oscillators. Therefore in the next section we analyze the relation between the 
frequency of oscillations and the synchronization in an array of harmonic os- 
cillators. The results of next section will be similar to the results of previous 
sections, however there will be two main differences. The first difference is that, 
unlike Lienard oscillators, where the magnitude of oscillations at synchroniza- 
tion is fixed, i.e., independent of initial conditions, with harmonic oscillators 
synchronization can occur at any magnitude depending on the initial condi- 
tions. The second difference is that the initial phases of harmonic oscillators do 
not play any role in determining whether synchronization will take place or not 
provided that the frequency of oscillations is large enough. Recall that this was 
not the case with Lienard oscillators. 



7 Synchronization of harmonic oscillators 

Consider the following array of coupled harmonic oscillators 

Qi = i^Pi (26a) 
Pi = -w^i + ^7ij(pj -pi) , i = l,...,m (26b) 

where, as before, w > is the frequency of oscillations and {7^^} is a connected 
interconnection. Functions 7^ are assumed to be locally Lipschitz. We let 
= [li Pi]'^ denote the state of zth oscillator. Array (p6| defines the below 
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system in R^™ 



e = h{^, io) (27) 



where ^ = [^f . . . ^m]"^ and what h is should be clear. Following the same 
procedure adopted for Lienard array, applying change of coordinates Xi{t) = 
e-^(")*6(0 yields 



E 



- sin ut 
cosujt 



7ij([— sinwi cosujt\(xj — Xi)) (28) 



which defines the below system in R^™ 

x = h'^'^{x,ujt) (29) 

for a; = [xj . . . a;^]^. Then, due to the periodicity of righthand side of we 
can talk about the average array 

~Vt) (30) 

where %j is as defined in ([7]). Array pO|) defines the below system in M^™ 

f, = h^'\^) (31) 

for r]^[ril ... ri^Y . 

In the remainder of the section we establish via sequence of three theorems 
the relation between frequency of oscillations and synchronization of coupled 
harmonic oscillators. In the first of those results (Theorem HI) we show that the 
oscillators of average array (|30p globally synchronize. Then from our first result 
we deduce (in Theorem [7]) that solutions Xi{-) of array (^5)) should eventually 
become arbitrarily close to each other, while initially being arbitrarily far from 
each other, provided that frequency of oscillations is arbitrarily large. Finally we 
claim (in Theorem El) that what is true for array ([28]) is also true for array p6p 
due to the norm-preserving nature of the transformation being used to transition 
between two arrays. 

Let us now define synchronization manifold A C K^™ as 

^ := {?7 e R^™ : rj, ^ rj^ for all i, j} 
to be used in the theorems to follow. 

Theorem 6 Consider system pip . Synchronization manifold A is globally 
asymptotically stable. 

Proof. By Lemma [2] we can write 



21 



for i = 1, . . . , m. Let G denote the graph of interconnection {7^ }. Graph Q is 
connected by assumption. Note that pij is continuous and zero at zero. Also, 
by Claim [21 if there is no edge of G from node i to node j then pij{s) = 0. If 
there is an edge from node i to node j then Pij{s) > for s > 0. 

Therefore hf~' is continuous; and vector hf~'{rj) always points to the (rela- 
tive) interior of the convex hull of the set {rii}U{rij : there is an edge of G from node i 
to node j}. These two conditions together with connectedness of G yield by [HI 
Cor. 3.9] that system ([9]) has the globally asymptotic state agreement property, 
see 10 Def. 3.4]. Another property of the system is invariance with respect to 
translations. That is, h}^'^{r] + C) = h^'^{ri) for ( E A. These properties let us 
write the following. 

(a) There exists a G /C such that |?7(i)|.A < a(|r/(0)|^) for all t >0. 

(b) For each r > and e > 0, there exists T > such that |77(0)|^ < r implies 
\vit)\A < e for ah t > T. 

Finally, (a) and (b) give us the result by [ini Prop. 1]. ■ 

Theorem 7 Consider system (|29p . Synchronization manifold A is semiglobally 
practically asymptotically stable. 



Proof. Consider array (|28)) . 
the relative distances Xj — Xi 
Let us define 



Observe that the righthand side depends only on 
This allows us to reduce the order of the system. 



X2 

2:3 



Xl 
Xl 

- Xl 



Note that y = h'{y, cot) for some h' : R2m-2 ^ M.>o -> ]R2m-2_ gj^^^g functions 
"fij are assumed to be locally Lipschitz, h^ is locally Lipschitz in y uniformly 
in t. Also, h'^ is periodic in time by p8|. Now consider array (l30l) . Again the 
righthand side depends only on the relative distances r/j — rji. Define 



m - m 
Vrn - m 



Then z — h^{z) where /i'' is the time average of h'' and locally Lipschitz both 
due to that is the time average of 7^ j . 

Theorem ini implies that the origin of i = h'^iz) is globally asymptotically sta- 
ble. Then [18l Thm. 2] tells us that the origin oi y = h'^{y, ujt) is semiglobally 
practically asymptotically stable. All there is left to complete the proof is to real- 
ize that semiglobal practical asymptotic stability of the origin oiy = h^{y, ujt) is 
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equivalent to semiglobal practical asymptotic stability of synchronization man- 
ifold A of system U 



Recall that system is obtained from system (^7]) by a time-varying 
change of coordinates that is a rotation in R^. Since rotation is a norm- 
preserving operation and set A is invariant under rotations, Theorem [7] yields 
the below result. 

Theorem 8 Consider harmonic oscillators (|27p . Synchronization manifold A 
is semiglohally practically asymptotically stable. 

8 Conclusion 

We have shown that nonlinearly-coupled Lienard-type oscillators (almost) syn- 
chronize provided that their phases initially lie in a semicircle and the frequency 
of oscillations is high enough. We have generated the same result for nonlinearly- 
coupled harmonic oscillators without any requirement on their initial phases. 
We have employed averaging techniques to establish our main theorems. 
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